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The T]' meson from lattice QCD. 
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Abstract 



We study the flavour singlet pseudoscalar mesons from first principles 
using lattice QCD. With Nj = 2 flavours of light quark, this is the so- 
called r]2 meson and we discuss the phenomenological status of this. Using 
maximally twisted-mass lattice QCD, we extract the mass of the r]2 me- 
son at two values of the lattice spacing for lighter quarks than previously 
discussed in the literature. We are able to estimate the mass value in the 
limit of light quarks with their physical masses. 



1 Introduction 



There is considerable interest in understanding hadronic decays involving rj and 
7]' in the final state. The phenomenological study of hadronic processes involving 
fiavour singlet pseudoscalar mesons makes assumptions about their composition. 
Here we address the issue of the nature of the r] and rj' from QCD directly, making 
use of lattice techniques. 

Lattice QCD directly provides a bridge between the underlying quark descrip- 
tion and the non-perturbative hadrons observed in experiment. The amplitudes 
to create a given meson from the vacuum with a particular operator made from 
quark fields are measurable, an example being the determination of /^r. It also 
allows a quantitative study of the disconnected quark contributions that arise in 
the fiavour singlet sector. The lattice approach provides other information such 
as that obtained by varying the number of quark fiavours and their masses. 

The disconnected diagram responsible for giving the fiavour-singlet pseudo- 
scalar meson a mass (in the chiral limit) is closely related to the fiuctuations in 
the topological charge. This link will be discussed shortly in this paper and in 
more detail in a later publication [1]. 

From the chiral perturbation theory description (for a review, see ref. [2]), one 
expects the mixing of rj and rj' to be most simply described in a quark model basis. 
In the fiavour singlet sector, for pseudoscalar mesons, we then have contributions 
to the mass squared matrix with quark model content (uu+dd) / \pl and ss (which 
we label as nn and ss respectively): 



Here m corresponds to the mass of the fiavour non-singlet eigenstate and is the 
contribution to the mass coming from connected fermion diagrams while x corre- 
sponds to the contribution from disconnected fermion diagrams. Thus rtinn is the 
pion mass. Because of mixing, as will be discussed, m^^ does not correspond to 
any specific meson, but its mass can be estimated assuming that the non-singlet 
pseudoscalar mass-squared is linear in the quark mass. So m^^ = 2m^^ — m^^, 
that is 2m|; — m^, leading to ra^s = 0.687 GeV. For a discussion from lattice 
results of small corrections to this assumption, see ref. [3]. Chiral perturbation 
theory also gives corrections to linearity coming from loop corrections. 

The mixing between the nn and ss fiavour singlet channels must produce the 
experimental ?7(548) and ?7'(958). One approximation used historically is that 
of SU(3) fiavour symmetry for which the two physical states will be a fiavour 
octet {nn — ^/2ss) and a fiavour singlet {\/2nn + ss). The t] is then identified as 
primarily the fiavour octet while the rj' is the fiavour singlet. 

Using as input m„„ and nigs and requiring that the output masses (m^ and 
m^/) are correctly reproduced, the three mixing parameters x cannot be fully 
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determined. It is usual to express the resulting one parameter freedom in terms 
of a mixing angle, here defined by 

?7 = rinn COS - riss sin , r/' = rinn sin (j) + r^ss cos (j) . (2) 

In our lattice study, we have two degenerate light quarks in the sea but no explicit 
strange quark. Since we do not consider partially quenched QCD here, we only 
study the top left corner of the mixing matrix. The flavour singlet pseudoscalar 
meson (called 772 in this case) will then have mass-squared m^^ + 2x„„. 

From the phenomenological analysis of the full mixing matrix, one can then 
estimate the mass of the 772 meson. One such phenomenological analysis, moti- 
vated by lattice input (basically the magnitude of x for strange quarks), gave [3] 
values Xnn = 0.292, Xns = 0.218, Xgs = 0.13 GeV^. This assignment corresponds 
to a mixing angle close to 45° and to a mass of the 772 meson of 0.776 GeV. 
This mass value is plausible since it is close to the (mass-squared weighted) av- 
erage mass of the 77 and 77'. Moreover, it is somewhat lighter than the t]' meson, 
which has an additional contribution to its mass from strange quark loops {xss)- 
However, here we plan to determine the 7/2 meson mass more directly, avoiding 
some of the assumptions made above. 

Note that the mass of the 772 meson comes from two components: rUnn which 
decreases with decreasing quark mass and Xnn which, in the phenomenological 
study, increases with decreasing quark mass. It is thus important to study the 772 
meson at light quark masses to explore this. For instance, with the parameters of 
the phenomenological model described above, the 7/2 meson made of light quarks 
of mass equal to the strange quark mass would have mass-squared m^g + 2xss 
giving mass 0.862 GeV. Thus we can expect a rather flat dependence of the 7/2 
mass on the mass of the underlying quarks (here a degenerate pair of valence 
quarks with sea quarks of the same nature). 

Lattice QCD is able to access fiavour-singlet states, but at a cost which has 
limited these explorations in practice. Since the correlator of a flavour singlet 
meson created at x and annihilated at y will have two contributions: a con- 
nected contribution where both quark and antiquark propagate from x to y and 
a disconnected contribution where there is a quark loop at x and another at 
y. This latter contribution needs to be measured at many values of x and y, 
so it is optimal to use a stochastic method to achieve this. One then evaluates 
the disconnected contributions (loops) at all x, and from combining them pair- 
wise, one can evaluate the disconnected contribution to the meson correlator. 
The stochastic method introduces noise, which can be reduced by appropriate 
variance reduction techniques. The goal is to reduce the noise to a level lower 
than the inherent noise coming from the underlying gauge-field variation. We 
shall show that the formalism known as twisted mass lattice QCD allows a very 
effective variance reduction method to be applied. 

We present details of our lattice methods and explain the stochastic method 
used to evaluate the flavour singlet pseudoscalar meson correlators (both con- 
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nected and disconnected). Because the signal from the disconnected correlator 
is relatively noisy, we discuss several strategies for reducing this error. From a 
combination of methods, we are able to present quite precise results, which enable 
us to discuss the continuum and chiral limit of the flavour singlet pseudoscalar 
mass. We compare with previous determinations and summarise. 



2 Twisted mass lattice QCD action 

In the gauge sector we employ the so-called tree-level Symanzik improved gauge 
action (tlSym) [1], viz. 



6„ ^ {1 - ReTr([;';;j} + 6, ^ {1 - ReTr(C/i;j,)} 

fl,U = l fJ.,U=l 



with the bare inverse gauge coupling /3, h\ = —1/12 and = 1 — 8bi. The 
fermionic action for two flavours of maximally twisted, mass degenerate quarks 
in the so called twisted basis [HI E] reads 

Stm = a'^^{x{x) [Dw[U]+mo + ifiq-f5T^]x{x)} , (3) 

X 

where rriQ is the untwisted bare quark mass, fig is the bare twisted quark mass, 
is the third Pauli matrix acting in flavour space and 

Dw[u] = l[i,{v, + v;)-av;v,] 

is the mass-less Wilson-Dirac operator. and V* are the forward and backward 
covariant difference operators, respectively. 

The mass term mo is tuned to maximal twist (as described in ref. [7j) at our 
lightest fig parameter. This guarantees 0{a) improvement [Ej. This was shown 
to work excellently in the quenched approximation [H [9], and there are good 
indications that this holds true also with dynamical fermions fTUl [TTj . More- 
over, this formalism has been found to give a very effective way to reach light 
quark masses [3, [121 [13] • As such, it provides an attractive route to evaluate the 
properties of the flavour- singlet pseudoscalar mesons. 

There is one complication, however, namely that twisted mass lattice formal- 
ism breaks flavour and parity symmetries at finite value of the lattice spacing 
a. These symmetries are restored in the continuum limit and the theory is well 
defined (so providing a valid regularisation) at finite lattice spacing. For our 
present purpose, this flavour-breaking causes the vr"*" and tt" mesons to have a 
mass splitting (of order a^). Moreover the 7r° has contributions from discon- 
nected diagrams [IH[7j. The analysis of the flavour- singlet pseudoscalar meson 
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Ensemble 




X T 


/5 


ajjLq 


1% 


ro/a 


Bi 




X 48 


3.9 


0.0040 


0.160856 


5.22(2) 


B3 








0.0085 








32^ 


X 64 


3.9 


0.0040 


0.160856 




Ci 


32^ 


X 64 


4.05 


0.003 


0.157010 


6.61(3) 


C2 








0.006 







Table 1: Summary of ensembles produced by the ETM collaboration used in this work. We 
give the lattice volume x T and the values of the inverse coupling (3, the twisted mass 
parameter a/iq, the hopping parameter k = (2aTOo + 8)""'^ and the Sommer parameter rg/a in 
the chiral limit from ref. [3 En]- The data sets cover 5000 equilibrated trajectories (10000 for 

is not significantly affected by this and it can be studied in the same way as was 
used previously in Wilson-based lattice formalisms [3, [T3]. In the twisted mass 
formalism, the 772 meson can mix in principle with the neutral Oq meson. This 
is an order a mixing which could induce an order contribution to the 772 mass 
observed at finite lattice spacing. We expect this effect to be negligible, both 
because the Oq meson is heavy and because we are working at maximal twist. We 
check this by varying the lattice spacing a. For a recent review see ref. |16j . 
For later convenience we introduce the following notation 

Du,d = Dw + mo ± 2/ig75 (4) 

for the fermion matrix of the two degenerate quark flavours (here labelled u and 
d) separately. 

The results presented in this paper are based on gauge conflgurations as pro- 
duced by the European Twisted Mass collaboration (ETMC). The details of the 
ensembles are described in ref. (TU]. In particular, we concentrated for this pa- 
per on the ensembles labelled Bi,B^,Bq and Ci and C2. We have compiled the 
details for those ensembles in table [H The 5-ensembles correspond to a value of 
the lattice spacing of about a ~ 0.09 fm and the C-ensembles to a ~ 0.07 fm. 
The spatial lattice size is of about L ~ 2.2 fm for all ensembles used here, apart 
from Bq, which has identical parameters to Bi but L ~ 2.7 fm. 

3 Neutral particles in twisted mass QCD 

In quenched or partially-quenched lattice QCD, the disconnected contribution to 
the flavour singlet meson does not combine properly with the connected contri- 
bution to give a physical state. To avoid this problem, it is mandatory to study 
flavour singlet states in full QCD - with sea quarks having the same properties 
as valence quarks. Then the spectrum of flavour singlet states is well deflned and 
can be extracted from the t-dependence of the full correlator. Here we focus on 
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the case where there are two degenerate hght quarks (called A'^^ = 2) which is a 
consistent theory in which to study the flavour singlet mesons. 

One promising way to study light quarks is using twisted mass QCD at maxi- 
mal twist. In the case of twisted mass fermions, the bilinear operator appropriate 
to create the 772 state is ip'^^ip which on transformation into the twisted basis (used 
in lattice evaluation) will become x^sXi where Ts acts in the (m, d) flavour space. 
This amounts to evaluating, in the lattice basis, the difference of the disconnected 
loop between u and d quarks. As already reported [TTl [TH], this enables a very 
efficient variance reduction technique to be used to evaluate the disconnected 
diagram relevant to the 772. 

The key observation is the following relation for the inverse and D^- 

l/D^ - l/Dd = -2z/i,(l/Z},)75(l/Z^„) . (5) 

Consider now the disconnected loop ^X(l/Du — l/Dd) where X is some 7- 
matrix (here J) and/or colour-matrix and the sum is over space. The conventional 
method involves solving 0^ = (1/-^")^^ with stochastic volume sources ^r- The 
desired result is then obtained from 

where the average is over noise samples (labelled r). However, the case mentioned 
above can be evaluated efficiently using the 'one-end-trick' [T9l [20] . Then the 
required disconnected loop is given by 

^X(l/D„ - = -2z/i,^(0*X750). . 

This has signal/noise which has a volume dependence V/ = 1 which is much 
more favourable than the conventional method with signal/noise 1/VV (here 
V = L^T). For ensemble Bi, we find the zero-momentum disconnected loop at a 
given time-value has a standard deviation of a = 18 for the variation with gauge 
configuration and time-slice (here called the intrinsic variation and obtained by 
extrapolating to an infinite number of stochastic samples) whereas the stochastic 
noise on this loop has a standard deviation of cr = 87 from 24 samples of volume 
source (conventional method) but only a = 7.5 from the above method (with 
12 samples). The relevant standard deviation for any analysis is the folding of 
the intrinsic variation [a = 18) with the stochastic error. So with a stochastic 
error of cr = 7.5 the net standard deviation is 19.5 whereas with the conventional 
method it would be 89. Thus there is a significant improvement. So 12 inversions 
give the disconnected correlator from all t to all t' with no significant increase in 
errors from the stochastic evaluation. 

We also use a non-local source/sink for the meson, constructed using "fuzzed" 
links of length 6a as described in refs. [2T1[T8]. which can be evaluated by replacing 
X by the corresponding fuzzed gauge links. Thus twisted mass lattice QCD allows 
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/5 


a/ig 


L/a 


No. 


am{ri2) 


xVdof 


3.90 


0.0040 


24 


888 


0.47(8) 


3/(24-6) 


3.90 


0.0040 


32 


184 


0.37(6) 


10/(24-6) 


3.90 


0.0085 


24 


245 


0.41(7) 


8/(24-6) 


4.05 


0.0030 


32 


198 


0.34(7) 


10/(24-6) 


4.05 


0.0060 


32 


188 


0.37(5) 


11/(24-6) 



Table 2: Results for the flavour singlet pseudoscalar meson from zero momentum correlators. 
The number of configurations analysed for neutral correlators (connected and disconnected) is 
shown (No.). The fits are to t-range 3-10 with 2 states and a 2 x 2 matrix of correlators (local 
and fuzzed at source and sink). The excited state mass is around am' = f .3 at /3 = 3.9. 

a very efficient way to evaluate flavour- singlet pseudoscalar disconnected loops, 
and combining them, the required disconnected correlator D{t). 

Evaluating the connected correlator C (t) for the neutral pseudoscalar meson 
as described elsewhere [18], then we have the full information C{t) — 2D{t) needed 
to explore the spectrum in the r]2 channel. We have a 2 x 2 matrix of such 
correlators available (local or non-local at sink/source). We could also explore 
pseudoscalar correlators obtained by creating the 772 state with V'7o75V' which 
on twisting becomes X7o75X- This does not allow our very efficient variance 
reduction to be applied, so the disconnected contributions are evaluated using 
a hopping parameter variance reduction [221 HE]- Even so, they are sufficiently 
noisy that we have been unable to use them to constrain the fits or to determine 
the r]2 decay constant. 

We measure the auto-correlation versus trajectory number for our correlators. 
The largest auto-correlation is found for the disconnected contribution. To ex- 
plore this more fully, we investigated the auto-correlation for the time-slice sum 
(i.e. zero momentum) of the quark loop at a given time versus trajectory, finding 
comparable auto-correlation times to that of the plaquette [10]. This is not un- 
expected, since the disconnected contribution is related to the topological charge 
density and thus both it and the plaquette are sensitive to the vacuum structure 
encoded in the gauge configurations. For further discussion see ref. [18]. In order 
to cope with this measured auto-correlation, we block the data into sufficiently 
large blocks (more than 80 trajectories) to remove any effect on our final error 
estimates. Note that the gauge configurations we use for our measurements are 
usually well separated in units of trajectories, such that even after blocking 80 
trajectories we are left with a sufficiently large number of (then independent) 
measurements. 

To explore the signal, we plot the effective mass versus t for ensemble Bi in 
fig. [H where we have used the variational basis from t-values 3 and 4 to optimise 
the ground state contribution. The zero momentum results (open symbols) show 
a plateau, although statistical errors are large, as will be discussed later. The 
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2 3 4 5 6 7 8 

t/a 

Figure 1: 772 effective mass (from a variational basis) versus t. Results are from ensemble Bi 
[L = 24) and {L = 32) (both a fig = 0.004). The results from momentum fc = 1 (in units 
of 27r/L) have been plotted as mass values assuming E'^ = + fc^. The results are slightly 
displaced for better readability. 



results for the two lattice spatial volumes available in this case agree well. 

In order to explore options to reduce the errors, we also evaluate the full 772 
correlator for momentum 1 (in units of 2tt /L). Since we evaluate the disconnected 
correlator for momentum 1 in each of the 3 spatial directions, we obtain an 
improved estimate of D (for example with relative error reduced to 11% at t = 10 
compared to error of 14% for zero-momentum). Since the symmetry classification 
of mesonic states is less sharp at non-zero momentum, the interpretation of the 
lattice spectrum is thus less straightforward. In practice, we find that the small 
reduction in error on D does not translate into a significant reduction in the error 
in determining the underlying mass of the 772 state. 

We make a fit to the combined 2x2 matrix of zero-momentum correlators 
for a t-range of 3-10 and with two states. The results from an uncorrelated fit to 
blocked data are presented in table O These fits are stable to varying the t-range. 

Before analysing these results further, we now discuss why the errors are so 
large for the 772 disconnected correlator, despite the fact that we measure all t 
and t', we use many gauge configurations and our stochastic errors are small. 
As described above, we have investigated the autocorrelation (versus trajectory) 
of the disconnected contribution and find no statistically significant evidence of 
any autocorrelation beyond 40 trajectories. We block our data into larger blocks 
(typically covering 80 trajectories) to avoid any increase in error from this source. 

The origin of the large error is that the signal for the disconnected part of 
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the correlator comes from only a small part of the total data sample. Here we 
concentrate on the zero- momentum method to illustrate this. We study Di{t' — t) 
where i labels the gauge configuration and we explore the distribution versus i 
at fixed \t' — t\. For instance (for ensemble Bi from 48 t- values for 888 gauge 
configurations) with \t' — t| = 10, 2.1% of the gauge configurations contribute 
26% of the observed signal (i.e. the sum over all gauge configurations). Thus the 
statistical impact of the data set is smaller than expected since parts of the data 
have big fluctuations (in a fermionic loop related to topological charge density). 

Another way to illustrate this error problem is that from the first 444 gauge 
configurations (covering 4440 trajectories) we find D/C = 0.96(29) at \t'—t\ = 16, 
while for the second 444 gauge configurations we find D/C = 0.47(23). Positivity 
requires D/C < 0.5. So the error estimate for certain quantities from even such 
a big ensemble as 5000 trajectories can be underestimated compared to having 
10000 trajectories, while other quantities do not show such a problem on the 
same set of gauge configurations. 

So even more configurations, than we have here, would be needed to get 
reliable and small errors in the case of disconnected contributions. This same 
conclusion has been obtained before, most strongly in a study involving staggered 
fermions |23j . 

4 Error reduction strategies 

Because the statistical error on the 772 mass is still large, we now discuss various 
strategies to reduce it while retaining the same set of gauge configurations. Some 
of these strategies do indeed reduce the statistical error, but at the expense of 
introducing a systematic error that has to be discussed. 

4.1 Excited state removal 

We now consider replacing the connected neutral pseudoscalar correlator C by 
just the ground state contribution. This has been considered previously [2l] . The 
basic idea is that the disconnected contributions {D) are only big for the lightest 
flavour- singlet pseudoscalar (the 772) and not for excited states. In other words, tt 
and ?72 are split but vr' and 773 are almost degenerate. This might happen because 
the topological charge fluctuations in the vacuum, which give the flavour singlet 
states a different mass, are more strongly coupled to the ground state 772 than to 
its excited states. Note that this is an assumption which needs to be checked by 
looking at the lattice results. 

In order to have no significant excited state contributions in the appropriate 
flavour- singlet correlator C — 2D, one should remove them from C. Then neither 
C nor D will have excited state contributions, and a one state fit to C — 2D should 
be possible down to quite low t-values, as indeed we shall find. The extension 
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/5 




L/a 


am{ri2) 




3.90 


0.0040 


24 


0.37(5) 


0.2/(27-3) 


3.90 


0.0040 


32 


0.30(4) 


2.4/(27-3) 


3.90 


0.0085 


24 


0.36(5) 


3.0/(27-3) 


4.05 


0.0030 


32 


0.26(5) 


1.8/(27-3) 


4.05 


0.0060 


32 


0.26(3) 


0.8/(27-3) 



Table 3: Results for the flavour singlet pseudoscalar mass from zero-momentum correlators. 
The connected correlator is taken from the fitted ground state component. The fits are to 
t-range 2-fO with 2 states and a 2 x 2 matrix of correlators (local and fuzzed at source and 
sink) . 

of this argument to twisted mass lattice QCD is not trivial, since there are also 
disconnected contributions to the neutral pion itself. Near the continuum limit, 
however, the argument goes through as before. 

In detail, we use fits of the form ccosh(— m(T/2 — t)) to the zero momentum 
data for the connected neutral correlator in the t-range 10-23 to determine the 
mass m and the coupling c. Note that this mass m corresponds to a neutral pion 
only if it is interpreted as a (partially-quenched or mixed-action) Osterwalder- 
Seiler state [H]; we shall call it the 'connected pion'. We then use our best fit 
parameters of m and c to construct the connected pion correlator, and we assign 
it with an error using the bootstrap method. This constructed connected pion 
correlator is finally used together with the disconnected contribution to build the 
full ?72 correlator, which we use even below t = 10. As shown in fig. [2], the effective 
mass is now essentially flat. Hence our result is consistent with dominance by 
a single 772 state down to surprisingly small t-values, where the data have small 
statistical errors. In contrast, fig. [T] shows the effective mass obtained without 
this excited state removal assumption, and it is less flat and has larger errors, 
although the plateau value is consistent. 

With the excited state removal assumption, we can then make one state fits 
to the resulting 2x2 matrix of correlators with t-range 2-10 and the results are 
reported in table O This shows that we achieve a significant reduction in the 
error on the 772 mass. This reduction comes at the cost of a possible systematic 
error coming from the assumption we have made. 

4.2 Point to point correlators 

We discussed above the possibility of using correlators at smaller t to reduce 
errors. Another approach, which we now discuss, is to reduce errors at all t- 
values by focussing on a different quantity. The method we describe has been 
used recently [25] in a study of the disconnected contributions to the rjc meson. 
It is mandatory to use an all-to-all method to estimate the disconnected loops 
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V V I' I' I ■'I^hL.-^hL-- 

LL I — ■ — I I 

LF I — e — , 
FF I — " — I 
POINT I — V — , 

2 4 6 8 10 

t/a 

Figure 2: Effective mass for zero-momentum 772 correlator for ensemble Bi {afiq = 0.004) 
taking as the connected contribution the ground state only. The horizontal line is the value 
given by a fit to the correlator data in the i-range 2-10. Also shown are results with point-to- 
point correlators with r = 2a, as described in the text. The results arc slightly displaced for 
improved readability. 

(i(x, t) at all X and t. The disconnected correlator of interest reads 

D(x-y,t-0 = Kx,t)rf(y,t')), 

assuming that (i(x, t) has no vacuum expectation value, as is the case for the 
flavour- singlet pseudoscalar meson because of the combined flavour-parity sym- 
metry of TMQCD. The approach we have considered above is to sum D over the 
spatial coordinates x and y. This has the advantage that the zero momentum 
correlator is studied which has the simplest theoretical interpretation: as a sum 
of exponentials for each state. Now D is actually peaked when |x — y| ^0 and 
is small at large |x — y|. Moreover we find that the absolute error is largely in- 
dependent of |x — y|, thus the relative error is smallest when |x — y| 0. We 
illustrate the dependence of D on r = |x — y| in fig. El The zero momentum 
contribution is a sum over this r-distribution with weight approximately and 
this enhances the larger r region which has larger relative statistical errors. Ba- 
sically the zero-momentum sum picks up (inherent gauge-time variation) noise 
by summing over large r where the signal is unimportant but the noise is still 
significant. 

For example, we find for ensemble Bi with \t — t'\ = 10 that the relative error 
on D when r = (x = y) is 50% of the relative error on D summed over both 
X and y. Thus the peak signal to noise is twice as big as that summed over 
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Figure 3: Disconnected contribution to the 772 correlator at time t and spatial separation 
r from point source and sink. Results are from ensemble Bi {afiq = 0.004) from 319 gauge 
configurations. The curves arc from the model fit described in the text. 



relative spatial separation. It might be thought that taking r 7^ would give 
an even smaller relative error, since 6 spatial directions are averaged (for spatial 
separations on axis). Although we find less error reduction than for r 7^ 0, 

we do make use of a further small error reduction and we select r = 2a (on axis). 

We now discuss how to extract the conventional physics (i.e. mass values) 
from correlators at fixed r (rather than summed over all r). Consider a correlator 
corresponding to a meson of mass m. In the region of t where the ground state 
contribution dominates, we expect the correlator to be dependent explicitly on 
the meson mass when point operators are used at source and sink. Indeed one 
way to extract this behaviour is by evaluating the 4-dimensional lattice Fourier 
transform: 

to model the data, where, to respect the periodic boundary conditions (for a 
meson), = 27ini/L and po = 27rno/T. In eq. El we take = 2sin(pj/2) 
which corresponds to the most local lattice realisation of the derivative in the 
Klein Gordon action. Other implementations of the derivative would shift the 
mass by corrections of order (ma)^ as discussed long ago [26j. We have checked 
that this expression correctly reproduces the r-dependence for the charged pion 
correlator in the t-region where the ground state dominates, using the mass value 
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previously determined from fitting the zero momentum correlator. Thus, though 
less familiar, one can use the expression of eq. Oto extract the mass value from the 
t-dependence at fixed r, in place of the usual expression c(e~™* + e^™''^^*^)/ (2m) 
used at zero momentum. 

To exploit this error reduction for the 772, we combine the disconnected con- 
tribution (with r = 2a as discussed above) with the neutral pion connected 
contribution also evaluated with r = 2a. This we evaluate using point sources. 
However, it is optimal, as discussed in the preceding section, to use only the 
ground state contribution to the connected neutral pion correlator, as we now 
discuss. 

For this ground-state connected component, we use eq. [6] with as input the 
mass m and coupling c determined from the zero-momentum fit to the neutral 
(connected) pion. The dependence of the resulting 772 correlator on t (at fixed 
r = 2a) can then be parametrised again by eq. El One way to illustrate this, 
as shown in fig. O is by plotting an 'effective mass' (here defined as the mass 
parameter that solves the dependence given by eq. [6] for two adjacent t- values). 
This shows that we again have a rather constant 'effective mass' even from low 
t-values, consistent with a description by one state only. Moreover, the errors 
on the 'effective mass' are smaller than from the fixed momentum methods used 
above. 

To determine the 772 mass by this method, we then fit the resulting combined 
?72 correlator to a two parameter expression (free mass and coupling) according 
to eq. El Since this approach is less familiar, we illustrate in fig. H] the fit versus 
t and also in fig. [3l the r-dependence of the disconnected contribution compared 
to the measured data. In detail, we use fits to the zero momentum data for the 
connected neutral correlator for the t-range 10-23 to determine the connected 
pion mass m and the coupling c, as above. From the bootstrap ensemble of 
values of m and c, we determine, at each t-value, the point-to-point connected 
pion correlator at fixed r = 2a with its associated error. This is then used together 
with the measured point-to-point disconnected contribution to construct the full 
ri2 correlator. The result is then fitted with one state using eq. El and a t-range 
2-10 to obtain the 772 mass value given in table HI 

Using different t-ranges gives the same mass value, within errors. Compared 
to using the zero-momentum approach described previously, the statistical error 
on the mass is much smaller. Here we have used local source and sink operators 
(since in this case the fixed r behaviour is given without extra parameters). What 
is especially helpful, however, is that the resultant data is well fitted by only one 
contribution for t > a. 

We have presented this new and powerful method for one ensemble (-Bi) but 
we also apply it (using r = 2a) to the other cases considered and the results are 
collected in table HI 

This point-to-point approach provides a useful reduction in the statistical 
error, though at the potential cost of introducing a source of systematic error 
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afiq 


L/a 


am{ri2) 




3.90 


0.0040 


24 


0.403(21) 


2.0/(9-2) 


3.90 


0.0040 


32 


0.365(19) 


0.1/(9-2) 


3.90 


0.0085 


24 


0.380(18) 


0.7/(9-2) 


4.05 


0.0030 


32 


0.302(24) 


0.0/(9-2) 


4.05 


0.0060 


32 


0.308(18) 


0.0/(9-2) 



Table 4: Results for the flavour singlet pseudoscalar mass using point-to-point correlators for 
the disconnected part and extracting the connected part from the fitted ground state compo- 
nent. 



in relating the quantity measured to that actually required. For local meson 
creation and destruction operators, the t-dependence at fixed r is given by the 
same parameters as the conventional zero-momentum case, up to possible non 
rotation-invariant contributions of order which we expect to be insignificant 
at large +r'^. For spatially smeared (or fuzzed) operators this 0(4) invariance 
is not present and we would need more parameters to describe the point-to- 
point correlators. We also find strong evidence of ground state dominance, which 
enables us to use the more precise data from smaller t-values. 



5 Summary of 772 mass 

In order to compare the error-reduced results from different lattice spacings, we 
plot them using the ro/a- values given in table[T]to create dimensionless quantities. 



We see in fig. 5(a) that the results are consistent with each other and are also 
consistent with the results presented above using the standard fixed-momentum 
analysis technique. This increases the statistical impact of our results. 

We also plot previous Nf = 2 results [271 CSl EH] that are in the quark mass 



and lattice spacing region we are exploring in fig. 5(b) Different groups using 
different lattice formulations obtain consistent results. Results from our two 
lattice spacings are consistent with each other and we are able to evaluate the 
flavour singlet mass closer to the continuum limit (and in an order improved 
formalism since we work at maximal twist) than hitherto. 

An extrapolation to the chiral limit of our results gives a value of the 772 mass 
of rom(?72) = 1.99(15) assuming a linear dependence. From our lattice results, we 
actually see no statistically significant evidence for any slope, and if we assume a 
constant behaviour, we would obtain rom(?72) = 2.00(5). Using tq = 0.454(7)fm 
obtained [7] from the ETMC evaluation of f^, we obtain an 772 mass of 0.865(65) 
GeV where the error reflects that from the linear extrapolation. This updates 
the phenomenological estimate of 0.776 GeV discussed above, which was in turn 
mainly motivated by lattice results at heavier quark masses than those we now 
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4 6 8 10 12 14 16 

t/a 

Figure 4: Singlet pscudoscalar correlator for ensemble Bi [jiq ~ 0.004) with local point sink 
and source. The correlator is illustrated for spatial separation r = 2a versus time t. The 
upper solid line is the connected contribution from the ground state neutral meson while the 
circles show the disconnected contribution. The squares are the net contribution to the r]2 and 
the behaviour from a one state fit to them is also shown. Note that the one state fit is not 
exponential, as discussed in the text. 



have available. We have additional systematic errors which are not fully under 
control: from the chiral extrapolation, the excited state removal assumption, 
the continuum extrapolation and any possible volume dependence. We do not 
see any sign that these systematic errors are required to our current statistical 
precision, but, as a conservative estimate, we consider these systematic errors to 
be comparable to the statistical error. So we quote an 772 mass of 0.865(65) (65) 
GeV. 

One of the main conclusions of our study is that we find strong evidence 
that the flavour singlet mass remains finite as the quark mass is reduced to the 
chiral limit. This has implications for the topological charge: it implies that the 
topological charge susceptibility must decrease as in the chiral limit [22] • We 
will discuss our results for topological charge more fully elsewhere [1]. 
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Figure 5: Comparison of values for the 772 mass versus quark; mass in units of tq w 0.45 fm. 
Values of the lattice spacing used are shown by listing ro/a. Some of the points have been 
displaced horizontally for legibility. The ETMC results shown in (a) are those from the 2-state 
fits to zero-momentum correlators (first method in text, open symbols) and from the error- 
reduction obtained by using point-to-point correlators and projecting to the ground state of 
the neutral pseudoscalar connected contribution (filled symbols) . In (b) we compare our results 
with reduced errors (filled symbols) to other Nf — 2 results from other collaborations [T51[?7l[^ 
(open symbols) with relatively light pions and small lattice spacing (a < 0.1 fm. ro/a > 4.5). 
The strange quark corresponds to (roTOjr)^ ~ 2.5. 
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